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INTRODUCTION 


It  is  becoming  increasing  evident  that  the  chemical  and  spectroscopic 
properties  of  molecules  containing  high  Z  atoms  can  be  calculated  using  relativistic 
effective  potentials  (REP).  The  REP  include  the  spin-orbit  and  are  functions  of 
the  total  orbital  angular  momentum,  j.  But  for  many  problems  in  determining  the 
electronic  structure  a  more  suitable  approach  is  to  decompose  the  REP  into  a 
non-relativistic  potential  (AREP),  in  which  the  spin-orbit  interaction  is  averaged 
away,  plus  the  spin-orbit  interaction  in  the  form  of  an  ab-initio  effective  spin- 
orbit  operator  (2,3,4).  Remarkably  enough  the  uranium  oxide  systems  described 
below  are  best  treated  by  determining  the  electronic  structure  using  the  relativistic 
AREP  operator  and  then  introducing  the  spin-orbit  interaction  with  the  AESOP 
operator  in  a  intermediate  coupling  calculation.  This  approach  limits  the  number 
of  states  interacting  through  the  spin-orbit  since  the  present  interest  is  in 
the  behavior  of  the  ground  state  and  low-lying  excited  states. 

The  goals  of  this  program  are  two-fold.  The  testing  and  development  of  the 
relativistic  effective  potentials  necessarily  forms  the  first  goal.  As  we  shall 
see  below  all  of  our  tests  to  date  have  been  successful.  This  allows  the  conconitant 


application  of  these  techniques  to  the  second  goal  which  is  to  determine  the 
spectra  of  the  uranium  oxide  postive  Ions.  A  report  of  the  spectra  of  U0+  was 
reported  at  the  DNA  meeting  on  the  Spectra  of  Uranium  Oxide  Species  held  at  the 
Las  Alamos  National  Laboratory  in  January  1983.  A  preliminary  note  was  also 
prepared  for  distribution  at  that  meeting  and  is  Included  as  Section  2  of  this 
report.  Both  the  Infrared  and  visible  spectra  of  U0+  was  predicted.  Strong 
fluorescence  pumping  in  the  red  part  of  the  visible  will  pump  excited  vibrational 


levels  of  the  ground  complex  of  states.  Under  low-pressure  excitation  conditions, 
then,  IR  vibrational  transitions  are  expected  around  925  cm  ^  as  well  as  a  weaker 
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group  of  electronic  transitions  which  can  range  from  about  200  cm”1  to  1300  cm-1. 

The  Intensities  of  the  electronic  transitions  will  depend  on  the  extent  of  the 
pumping  of  the  ground  electronic  states. 

Section  3  will  summarize  the  preliminary  assignment  of  the  FeO  spectra. 

The  assignments  are  still  tentative  since  quantitative  calculation  of  the  differential 
correlation  energy  is  still  in  progress.  The  present  result  do  account  for  all 
of  the  observed  levels  below  1.5eV. 

Section  4  summarizes  the  calculation  of  the  spin-orbit  interaction  as  a 
function  of  the  internuclear  distance  for  the  ground  state,  X  n,  of  OH.  This 
system  was  chosen  as  a  test  of  the  AESOP  since  very  accurate  experimental  and 
theoretical  data  has  just  become  available.  In  addition,  this  system  is  a 
prototype  for  the  spin  uncoupling  that  accompanies  dissociation  into  two  open- 
shell  fragments.  The  present  calculation  accurately  reproduces  the  shape  of 
both  the  experimental  and  theoretical  all -electron  spin-orbit  curves.  The  present 
calculation  is  a  very  simple  one,  though,  compared  to  the  all-electron  result 
and  is  readily  extended  to  heavier  systems  while  the  all -electron  calculation  can 
only  be  extended  at  great  expense. 

In  the  semi-annual  report  we  noted  that  several  of  the  test  calculations  on 
the  effective  potentials  have  been  accepated  for  publication  including  references 
3  and  4  on  the  AESOP  operator  and  a  paper  entitled,  "Electronic  Structure  of  UH, 

UF,  and  Their  Ions,"  which  has  been  accepted  by  the  Journal  of  Computational 
Chemistry  [5].  We  have  also  completed  a  manuscript  on  "Finite-Field  SCF  Calculations 
of  the  Dipole  Polarizabilities  of  Heavy  Atoms  Using  Relativistic  Effective  Potentials," 
which  has  been  submitted  to  the  Journal  of  Physics,  London  [6]. 
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Relativistic  effective  potentials  are  ised  to  calculate  the  electronic 
structure  and  spectroscopic  properties  of  U0H .  The  lowest  energy  states  are 
very  ionic  and  the  molecular  orbitals  substantially  localized  so  that  the 
molecule  is  described  by  the  ionic  fragments,  U+3(f3,^I)  and  0~2(^S).  All  of 
the  quartet  states  from  £“  to  I,  that  arise  from  the  coupling  of  these  ionic 
fragments,  are  perturbatively  mixed  using  an  effective  operator  for  the  spin-  ■ 
orbit.  The  R  and  «e  of  the  ground  f;=  9/2  state  have  been  determined  to  be 
3.48  bohr  and  925  enf^ .  The  vibrational  and  electronic  states  are  interleaved 
with  the  lowest  excited  state,  fi=7/2,  at  1315  cm“^ . 

The  excitation  energies  of  the  excited  states  of  U0+  were  calculated 
using  a  restricted  valence  configuration  interaction.  Strong  radiative 
transitions  are  predicted  in  the  red  part  of  the  visible.  These  transitions 
are  predominately  atomic-like  f  to  d. 
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INTRODUCTION 

The  study  of  the  positive  ions  of  uranium  oxide  dates  from  the 
observation  of  U0+  and  U02+  from  the  interaction  of  U  +  02  in  a  molecular 
beam  experiment  [1,2].  Although  the  infrared  spectra  of  matrix  isolated 
neutral  uranium  oxide  species  from  U0.to  U03  [3,4,5]  and  the  spectra  of  U02+ 
have  been  assigned  [6],  there  is  no  reported  spectrum  of  U0+.  In  order  to  assist 
in  the  search  for  the  IR  and  visible  spectra  of  this  ion,  we  have  calculated  the 
electronic  structure  and  energies  of  the  lowest  4I  and  4H  states  of  the  ion  and 
selected  excited  states  which  can  couple  radiatively  with  the  4I  ground  state. 

The  use  of  relativistic  effective  potentials  (REP)  makes  these 
calculations  possible.  Recent  work  on  low-Z  molecules  with  non-relativistic 
effective  potentials  has  shown  that  the  energy  curves  compare  well  with 
all-electron  ab  initio  curves  [7,8].  Lee,  et  al  [9]  have  described  the 
method  for  calculating  the  j -dependent  REP  based  on  pseudospinors.  The 

6Si/2,  f2 *  6p3^2»  ^^2’  **^7/2*  ^^3/2*  6dg^2  were  obtained 

from  a  Dirac-Fock  solution  [10]  of  a  degeneracy  weighted  energy  average  of 

2  6  3  +? 

the  (core)6s  6p  5f  6d  open-shell  configuration  of  U  and  the  6g  spinors 

2  6 

were  obtained  from  the  (core)6s  6p  6g  ionic  configurations.  The  spin- 
orbit  interaction  is  very  large  for  uranium,  but  the  large  components  of  the 
6d  and  5f  spinors  with  different  j  values  are  not  substantially  different. 

The  positions  of  the  peaks  of  the  outermost  maxima  of  the  5fg^2  and  5f^2 
orbitals  differ  by  0.020  bohr  and  the  6d3^2  and  6dg^2  differ  by  0.092  bohr.  The 
calculation  of  the  electronic  structures  of  the  ground  states  of  these  metal  oxides 
are  more  readily  done  initially  with  neglect  of  the  spin-orbit  interaction  using 
A-dependent  relativistic  effective  potentials,  and  subsequently  introducing  the 


spin-orbit  Interaction  through  configuration  mixing. 

The  spin-orbit  components,  or  j-dependence,  of  the  REP  is  eliminated  by 
taking  a  weighted  average  [11]  of  the  j=£+l/2  and  j=£-l/2  potentials: 


yAREP 


U 


REP 

£,£-1/2 


+ 


£+1  „  REP 
2£+l  u£,£+ 1/2 


The  average  potential  has  been  designated  AREP.  Using  the  ALCHEMY 

system  of  codes  [12],  modified  to  use  the  AREP,  we  performed  preliminary 

self-consistent-field  (SCF)  calculations  which  determined  that  the  electronic 

structure  of  the  lowest  energy  states  is  very  ionic,  and  the  molecular  orbitals 

are  substantially  localized.  All  of  the  lowest  energy  states  arise  from 

the  maximum  projection  of  the  angular  momentum  of  the  ionic  fragments,  U+3(f3,4I) 
2  1 

0"  (  S),  on  the  molecular  axis.  The  lowest  energy  state  at  the  SCF  level  is 

4 

the  I.  All  of  the  SCF  quartet  states  from  £  to  I  are  found  to  be  similar 
in  electronic  structure  with  the  f  orbitals  being  highly  localized.  Radiative 

4 

transitions  from  the  I  state  will  be  assumed  to  be  representative  of  the 
complex.  An  accurate  calculation  including  correlation  is  very  difficult  because 
important  correlating  configurations  are  of  different  ionicity  or  orbital 
occupancy  from  the  SCF  configurations,  and  differential  correlation  effects  for 
the  ions  embedded  in  the  molecule  complicate  the  calculation.  A  qualitative 
description  of  the  electronic  structure  and  spectroscopy  is  still  possible 
because  the  U  3(f  »  I)  state  is  energetically  separate  from  the  other  atomic 

3 

ion  states.  The  f  doublets,  for  example,  are  nearly  2  eV  higher  in  energy. 

The  spin-orbit  mixing  of  the  lowest  energy  4E~,  4n,  4A,  4$,  4r,  4H,  and  4I 
states  has  been  studied  perturbatively  using  the  AESOP  operator  [13].  The  lowest 
fl  states  are  still  dominated  by  a  single  A  function.  The  mixing  coefficient  of 

4 

the  I  function  in  the  fia9/2  ground  state,  for  example,  was  calculated  to  be 

4 

0.94.  Transitions  from  the  dominant  I  component  will  then  provide  a  zero  order 
description  of  absorption  from  the  ground  state.  Small  MC-SCF  calculations  were 
done  to  show  that  the  SCF  electronic  description  was  adequate  for  a  qualitative 
understanding  of  the  ground  state. 
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The  excited  state  calculations  were  done  by  restriced  valence  configuration 

interaction  (VCI)  using  several  sets  of  localized  molecular  orbitals.  Since  we 

are  only  Interested  In  those  states  which  are  dipole  radiatively  coupled  to  the 
4 

I  ground  state »  we  will  describe  in  this  note  the  electronic  structure  of 
the  low-lying  4I,4H,  and  4K  states.  These  data  will  be  used  to  qualitatively 
describe  the  fluorescence  pumping  of  cold  U0+.  These  calculations  are  presented 
to  assist  in  the  planning  and  interpretation  of  experiments  designed  to  detect 
U0+  *  and  to  present  a  preliminary  analysis  of  the  electronic  structure  of  such 
heavy  metal  oxides. 

DESCRIPTION  OF  THE  CALCULATIONS 

A  double  zeta  Slater  function  atomic  basis  was  energy  optimized  using 
the  AREP  on  the  U+2  (6s26p66d5f3,5L),  U+1(6s26p67s6d5f3,6L),  and  U+1  (6s26p67p6d5f3 
^M)  states.  SCF  calculations  on  the  4I,  4H,  and  6K  molecular  states  provided  a 
set  of  molecular  orbitals  primarily  localized  on  the  oxygen  pa  and  pw  orbitals, 
and  the  uranium  6s,6p,fa,f7r,f6,f4>,dir,d6,  and  7sa  orbitals.  In  order  to 
determine  the  spin-orbit  coupling  among  the  lowest  energy  4jf,  4n,  4A,  4$,  4r, 

4H,  and  4I  states  that  arise  from  all  f3  couplings,  it  was  necessary  to 

calculate  these  states  with  a  common  orthogonal  set.  Using  the  SCF  molecular 

4  4 

orbitals  from  the  I  and  H  states,  configuration  interaction  calculations 
were  carried  out  using  a  base  set  of  configurations  for  each  state  which  arise 
from  three-fold  occupancy  of  the  fa,fn,f6»  or  f<j>  orbitals  and  adding  all 
configurations  arising  from  a  single  excitation  from  all  valence  orbitals  in 
the  base  configurations  into  the  virtual  space.  The  spin-orbit  interaction 
was  calculated  using  the  AESOP  operator.  The  energies  of  these  $2  coupled  states 
are  given  in  Table  1  at  one  distance  (3.625  a.u.)  relative  to  the  energy  of 
the  unperturbed  4I  state. 
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The  lowest  energy  ft  states  are  still  dominated  by  one  A  state.  For 
example,  the  ft=9/2  ground  state  is  represented  by 

4Ig/2  =  0.94  4I  -  0.33  4H  -  0.10  4r  +  0.02  4<f> 

4 

and  the  H first  excited  state  is  given  by 

4h7/2  =  0,83  4”  "  0,48  4r  "  °-25'4^  +  °-08  4A- 

4  4  - 

The  energy  splitting  of  the  I  through  Z  states  is  about  0.3eV,  but  after 
considering  the  spin-orbit  coupling,  the  complex  of  ft  states  is  spread 

3 

over  leV.  Since  spin-orbit  coupling  to  states  outside  the  quartet  f  groups 
was  neglected,  only  the  description  of  the  lowest  energy  states  has  any 
validity. 

A  Morse  fit  to  the  4I  SCF  curve  yields  an  equilibrium  internuclear  distance 
of  3.46  bohr  and  a  vibrational  constant  of  949  cnf^  for  U*60+.  The  4Ig/2  9round 
states  energy  curve  is  found  to  have  the  same  Re  and  a  vibrational  constant  that 
is  only  10  cm”^  smaller.  A  small  multi-configuration  SCF  (MCSCF)  calculation 
for  the  4I  state  that  includes  the  dominant  back  transfer  correlation  yields 
an  R.  of  3.48  bohr  and  an  u>  of  935  cm"^.  Assuming  the  spin-orbit  correction 

G  ® 

to  the  4I  Cl  function  calculated  above  applies  to  the  MC-SCF  result,  the 
vibrational  constant  for  the  spin-orbit  perturbed  MC-SCF  energy  curve  is  935  cm  ^ • 

4 

A  restricted  valence  configuration  interaction  was  performed  on  the  I, 

4  4 

H,  and  K  states  to  provide  an  estimate  of  the  absorption  spectrum.  The  U6s, 

U6p,  and  02s  orbitals  were  retained  fully  occupied  and  the  fn,dn,  and  d6  orbitals 
were  obtained  In  the  6K  SCF  solution.  This  calculation  does  not  allow  for  orbital 

4 

relaxation,  so  the  energy  curves  are  not  optimal.  The  Re  calculated  for  the  I 
state  is  about  3.6  a.u.  in  this  calculation.  This  value  is  about  0.1  a.u.  larger 
than  that  obtained  in  the  MC-SCF  calculation  and  is  attributed  to  the  use  of 
orbitals  that  are  intended  to  describe  both  the  ground  and  excited  states  and  are 
not  optimal  for  either. 
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The  excitation  wavelengths  in  the  visible  region  from  the  ^1  ground  state 

are  shown  in  Figure  1.  All  the  curves  are  shifted  so  that  the  ^1  agrees  with 

the  value  calculated  for  the  MC-'SCF  curve  while  retaining  their  relative  positions. 

The  VCI  total  correlation  energy  near  Rg  is  about  leV  for  both  the  ^1  and  ground 

states  and  results  from  two  types  of  configuration  Interaction.  The  first  is  the 

double  excitation  between  the  slightly  bonding  (aU7s  +  02p)4a  valence  orbital,  which 

is  doubly  occupied  in  the  SCF  wavefunction  and  the  slightly  anti-bonding  (U7s-a02p)5a 

correlating  orbital  (a«l).  The  second  type  of  correlation  mixing  results  from  back- 

bonding  from  the  (02p)2ir  valence  orbital  into  tt  correlating  orbitals.  The 

4 

dominant  correlating  orbital  Is  primarily  U6dTr  in  character  for  I  and  a  mixture 

of  U5fir  and  U6dir  for  ^H.  The  same  correlating  configurations  were  included 

4 

in  the  MC-SCF  calculation  of  the  I  ground  state. 

4  4 

Transition  moments  from  the  ground  I  state  to  excited  I  states  and 
4  4 

all  the  H  and  K  states  are  given  in  Table  2  for  R=3.625  bohr.  The  four 
strongest  transitions  which  range  from  0.75D  to  over  ID  are  all  atomic  f  to  d 
transitions.  These  electronic  transitions  range  from  the  far  red  around  750  nm 

4 

to  the  green  at  455  nm.  The  electronic  transition  between  the  ground  I  and 
4H  states  has  a  relatively  weak  transition  moment  of  0.07  D. 

DISCUSSION 

The  ground  state  electronic  structure  and  spectroscopic  constants  are 
based  on  SCF  and  MC-SCF  calculations  using  the  AREP  with  inclusion  of  the 
spin-orbit  coupling  as  a  perturbation.  The  MC-SCF  calculation  shows  little 
change  in  the  orbitals  from  the  SCF  results  with  only  a  modest  shift  in 
the  dipole  moment  (M).2  a.u.)  as  a  result  of  the  it  back-bonding.  However, 
at  least  two  major  correlation  effects  have  been  neglected  because  of  the 
inability  to  handle  the  calculation  at  this  time.  The  differential  atomic- 
like  correlation  between  states  of  different  ionicity  and  d  or  f  occupancy 
and  charge  transfer  plus  polarization  configurations  must  be  considered  at 


some  future  date.  As  a  result  it  is  impossible  to  gauge  the  accuracy  of 
the  calculations  without  recourse  to  comparison  with  experiment. 

An  estimate  of  the  ground  state  vibrational  frequency  can  be  obtained 

by  comparing  the  spectroscopic  constants  between  the  SCF  curves  for  the 
4  +  5 

I  state  of  UO  and  the  I  state  of  UO.  The  VCI  curve  is  not  sufficiently 
reliable  because  the  correlation  is  not  optimized  as  a  function  of  the 
distance.  The  best  Morse  fits  to  the  SCF  curves  yield  an  equilibrium 
internuclear  distance  of  3.46  bohr  and  a  vibrational  constant  of  949  cm”1 
for  Ul60+  compared  to  3.56  bohr  and  863  cm-1  for  Ul60.  A  band  at  about 
820  cm"1  has  been  assigned  to  U160  in  the  Ar  matrix-isolated  spectrum  [5], 
but  this  value  may  be  appreciably  red-shifted  from  the  gas  phase  value.  For 
the  Kr  matrix  the  820  cm"1  band  is  split  and  red-shifted  to  819  and  815  cm"1. 
These  shifts  are  smaller  than  analogous  shifts  observed  in  the  lithium 
halides  [14].  Linevsky  [14a]  has  shown  that  the  dominant  interaction  causing 
the  shift  is  the  dipole  induced-dipole  coupling  of  nearest  neighbors  with 
the  difference  in  the  square  of  the  dipole  moments  in  the  v=0  and  1  levels 
determining  the  magnitude  of  the  shift.  Both  UO  and  the  alkali  halides  have 
large  dipole  moment  derivatives  which  we  find  are  of  comparable  magnitudes. 

But  the  calculated  value  of  the  anharmonic  constant  for  UO  is  about  2  cm"1 
or  about  half  the  value  for  LiCS,  [15].  The  variation  of  the  dipole  moment  as 
a  function  of  the  vibrational  state  in  UO  is  then  much  less  than  that  in  the 
alkali  halides.  However,  it  must  be  recognized  that  the  u>exe  is  obtained 
from  a  Morse  fit  to  an  energy  curve  that  is  not  represented  well  by  a  Morse 
curve.  The  UO  MC-SCF  calculations  yield  an  of  3.57  bohr  and  an  we  about 
845  cm"1  or  about  3%  larger  than  the  matrix  isolation  observed  frequency.  The 
spin-orbit  coupling  will  be  more  complicated  in  UO  than  U0+  but  the  additional 
states  that  mix  would  tend  to  lower  the  frequency.  We  thus  feel  that  the 
accuracy  is  within  the  3%  range  for  the  U0+  u>e  of  925  cm"1  determined  from  the 
spin-orbit  perturbed  MC-SCF  energy  curve. 


The  vibrational  infrared  transition  for  the  state  is  very  strong 

with  an  oscillator  strength  of  about  5.10"5  for  the  v=0  to  v=l  transition.  This 

compares  with  an  electronic  oscillator  strength  of  about  2.2.10"5  for  the 

4  4 

Ig/2  “  H7/2  infrared  transition.  The  vibrational  and  electronic  infrared 

transitions  will  interleave.  The  present  calculation  predicts,  for  example, 

the  lowest  energy  transition  to  be  a  vibrational  transition  at  about  923  cm"1 

followed  by  an  electronic  transition  at  about  1315  cm"1. 

The  visible  electronic  transitions  will  be  split  by  the  spin-orbit 

interaction  into  a  range  over  1  eV  wide  roughly  centered  around  the  transitions 

given  in  Fig.  1.  The  spin-orbit  interaction  for  the  excited  states  is  too 

complicated  to  be  treated  by  the  present  methods,  but  the  prediction  of 

strong  radiative  transitions  in  the  visible  will  not  be  affected.  The  quartet 

states  that  couple  to  the  lowest  4E,  4r,  and  4H  states  are  also 

2 

predominately  f  d  systems  and  have  the  same  transition  moment  behavior  as  the 
4  4  4 

H,  I,  and  K  excited  states.  Calculations  on  the  sextet  states  indicate  a  dominant 

+2  -1 

configuration  with  ionicity  U  0  '  that  will  not  couple  strongly  to  the  radiatively 

+3  -2  3 

important  quartets  with  ionicity  U  0  .  The  f  doublets  are  in  this  energy  region 

but  again  the  spin-orbit  coupling  is  weak  because  the  f  and  d  orbitals,  though 
polarized,  do  retain  considerable  atomic  character. 

Strong  radiative  transitions  are  expected  in  the  red  part  of  the  visible. 

The  fluorescence  will  pump  the  excited  vibrational  levels  of  the  lower 
complex  of  states  (which  have  similar  vibrational  structure)  because  the 
excited  states  are  shifted  in  frequency  and  Rg.  Under  low-pressure  excitation 
conditions  the  strong  visible  and  IR  transitions  at  about  925  cm”1  would  be 
expected.  Detailed  modeling  is  not  possible  because  of  the  true  complexity 
of  this  system. 
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TABLE  II 

4 

Transition  Moments  from  the  X  I  State  to  the  Excited  I,  H,  and 

Transition  Moment  (e  -  bohr) 

Roots 

State  1  2  3 

4I  0.31  0.092 

4H  0.028  0.31  0.61 

4K 


States 


0.30 


Absorption  Spectra  of  UO 


•H 


3.  Electronic  Structure  of  FeO 

A  survey  of  the  FeO  calculations  was  given  In  the  semi-annual  report.  The 
REP  was  shown  to  yield  energy  curves  comparable  to  all -electron  values  at  the 
SCF  level  where  comparisions  can  be  made.  The  configuration  interaction  calculations 
included  all  the  valence  bonding  correlation  energy  but  this  is  insufficient 
to  determine  accurate  energy  differences  because  of  differential  correlation  energy 
effects.  The  valence  states  can  have  different  lonlcity  or  a  different  number  of 
occupied  d  electrons  and  the  "atomic"  correlation  of  the  dominant  configurations 
is  then  required.  Such  calculations  are  presently  underway  but  we  will  show  that 
assignments  of  the  spectra  and  an  ordering  of  the  states  can  still  be  made. 

Interplay  between  the  ab  Initio  results  and  experimental  data  permits  the  precise 
description  of  the  electronic  structure  of  the  ground  and  excited  states  of 
FeO. 

5  5  +  5  5  7  + 

The  lowest  energy  X  A,  £  ,  $,  n,  and  £  states  will  be  analyzed  in  this 

7  + 

order.  Experimental  data  exsits  only  for  the  first  four  states.  The  £  state 
provides  a  good  example  of  the  obstacles  that  differential  correlation  errors 
raise  in  the  path  of  an  accurate  assignment.  This  state  is  actually  calculated 
to  be  lowest  In  energy  at  the  valence  configuration  interaction  (VCI)  level  but 
it  is  evident  from  the  orbital  populations  that  the  'atomic'  differential  correlation 
will  lead  to  this  state  being  over  leV  above  the  ground  state. 

A.  The  5A  State 

The  ground  state  of  FeO  is  now  known  to  be  the  5  A  from  unequivacal  experimental 
data  [1,2].  Earlier  calculations  had  also  suggested  a  5A  ground  state  [3,4]. 

5 

Cheung  et  al.  report  the  spectroscopic  parameters  of  the  A  state  and,  most 
Important  from  our  point  of  view,  they  determined  the  spin-orbit  splitting. 


ml 
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Because  the  differential  correlation  errors  result  both  in  errors  in 
calculated  excitation  energies  and  spectroscopic  constants,  we  must  look,  at  this 
stage  of  the  calculation,  for  other  properties  to  Identify  the  calculated  states 
with  the  observed  ones.  Cheung  et  al.  report  the  spin-orbit  coupling  constant, 
A=-94.9  cm”1.  A  complete  active  space  multi -configuration  self-consistent-field 
(CAS-MC-SCF)  calculation  was  performed  for  the  ^A  state.  The  02s  and  Fe3s  and 
3P  atomic-like  orbitals  were  kept  doubly  occupied  and  the  three  sigma,  two  pi, 
and  one  delta  orbitals,  that  span  the  02p  and  Fe4s  and  3d  orbitals,  were 

c 

populated  with  12  electrons  in  all  possible  ways  that  yield  A.  This  calculation 
yields  a  variationally  relaxed  valence  configuration  interaction  wave  function 
and  energy  for  the  ground  state  of  FeO.  Using  the  CAS-MC-SCF  wave  function  and 
the  AESOF  operators  for  Fe  and  0,  we  obtain  an  A=-85  cm"1  for  R=3.2  bohr.  This 
includes  a  Breit  correction  of  0.895  for  the  Fe  3d  spin-orbit  operator.  The 
error  of  about  10  cm"1  derives  primarily  from  the  neglect  of  the  residual  spin- 
other-orbit  valence  interactions  among  the  open-shell  d  orbitals.  These  interactions 
arise  from  the  Breit  interaction  and  are  not  Included  in  the  Dirac-Fock  atomic 
calculations  from  which  the  effective  potentials  are  derived. 

The  agreement  between  experimental  and  calculated  spin-orbit  constants  of 
around  10%  supports  the  assignment  of  the  5A  as  the  ground  state.  The  electronic 
structure  of  the  state  can  also  be  described  in  terms  of  the  dominant  configuration 
of  the  Cl  or  the  populations  of  the  natural  orbitals  that  represent  the  Cl  density 
matrix.  A  more  complete  description  of  these  functions  will  be  given  in  a 

2  2 

subsequent  report  but  the  dominant  configuration  in  this  case  is  (core)  la  2o  in 
4  2  2  4  2  3 

3o  4o  211  5o3n  16  ,  where  the  dominant  atomic  orbital  contribution  in  the 
molecular  orbitals  are: 

la  -  Fe3s,  2o-Fe3p,  UI-Fe3p,  3a-02s, 

4a-(Fe3d  +  02p),  2II-02p,  5a  -  Fe4s,  3n  -  Fe3d,  16  -  Fe3d. 
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The  most  significant  correlating  configurations  arise  from  excitation  from  the 
bonding  4o  orbital  to  the  anti-bonding  60.  Near  Rg  the  population  in  the 
60  orbital  Is  nearly  0.5. 

The  MC-SCF  wavefuntion  is  remarkably  different  from  the  SCF  wavefuntion. 

The  4o  orbital  in  the  SCF  function  is  almost  purely  02p  and  the  5a  is  composed 
primarily  of  Fe3d.  The  SCF  function  has  the  ionic  character  Fe+20~2  but  the 
back-transfer  correlation  allows  the  formation  of  a  partial  sigma  bond.  At  the 
relatively  short  equilibrium  Internuclear  distance,  the  stronger  bond  Is 
formed  with  the  Fe3d  atomic  function.  As  a  result  the  open-shell  5a  orbital 
transforms  from  Fe3d  in  the  SCF  function  to  Fe4s  in  the  MC-SCF  function.  This 
result  points  to  the  extreme  caution  required  in  accepting  the  result  of 
configuration  interaction  calculations  based  on  SCF  molecular  orbitals,  e.g., 
the  Cl  reported  for  MnO  [5]. 

B.  The  5I+  State 

Engelklng  and  Lineberger  [2]  have  found  an  excited  state  in  the 
photodetachment  spectrum  of  Fe0“.  Using  the  theoretical  calculations  for  both 

-  5  + 

FeO  and  FeO  ,  they  assigned  the  state  as  a  E  .  Our  calculations  support  the 

assignment.  A  low-lying  5E+  state  is  calculated  with  the  dominant  configuration 

(core)  lo22a2in^3a24a22II^5a2in2162,  but  the  molecular  orbitals  do  not  have  the 

5 

same  atomic  character  as  they  do  in  the  A  state.  The  natural  orbitals  for  the 

5  +  5 

MC-VCI  solution  show  that  the  E  state  has  far  more  s  character  than  the  A  but 

+2  -2 

both  states  are  predominately  Fe  0  .  Because  of  the  polarized  character  of 
the  s  containing  orbital,  the  correlation  energy  is  likely  to  be  less  for  the 
VCI  description  of  the  ^Z+  state  than  for  the  ^A.  The  MC-VCI  calculations  finds 
these  two  states  within  O.leV  of  each  other.  On  the  basis  of  the  differental 
correlation  energy,  the  ®E+  state  is  a  good  assignment  for  the  excited  state  at 
0.5eV. 
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The  ground  state  of  FeO"  is  calculated  to  be  a  ^A.  Both  the  ground  ^A 
and  the  excited  Z  states  can  be  obtained  from  this  anion  by  removing  a  single 
electron  from  an  orbital  in  the  dominant  configuration  of  this  state.  The  orbitals 
of  the  anion  are  substantially  different  In  their  atomic  structure  from  the 
neutral  states  but  the  overlap  Is  still  large  enough  to  rationalize  the  spectra 
observed  by  Engel  king  and  Lineberger. 

C.  The  5n  and  5$  States 

These  states  are  observed  by  Merer  et  al.  to  form  a  complex  of  levels  that 
lie  within  250  cm  "*  of  one  another.  The  closeness  of  these  levels  is  explained 
by  the  fact  that  the  dominant  configuration  is  the  same  for  both  states ,  (core) 
10^20^111^30^40^211^50111^16^60.  The  ionicity  of  these  states  is  predominately 
Fe+0~  with  the  0"  oriented  n  and  the  Fe+  in  a  predominately  sd®  ®D6  coupling. 

The  6o  is  dominated  by  a  very  polarized  Fe4s  orbital.  The  0"n  and  Fe+6  'atoms' 
are  weakly  coupled  to  produce  the  n  and  4>  states. 

The  spin-orbit  interaction  will  be  used  to  confirm  the  electronic  structure 
description.  Merer  has  observed  the  spin-orbit  splittings  between  the  ft  states 
of  the  n  and  $  states.  The  data  shows  significant  perturbations  in  both  the  n 
and  <i>  states  which  is  not  quantitatively  understood.  Estimates  of  the  spin-orjit 
coupling  constant  can  be  made  from  the  data,  though,  for  the  purpose  of  comparing 
with  the  calculated  vaules.  For  the  4>  and  n  states  the  experimental  coupling 

_  -1  c  c 

constants.  A,  are  47  and  200  cm  ,  respectively,  for  the  $  and  n  states.  The 
calculated  constants  are  41.9  and  195  cm"'*,  respectively,  for  the  5<i>  and  5n  states. 
The  comparison  is  good  enough  to  conclude  that  the  overall  electronic  description 
of  these  states  is  understood.  The  energy  curves  of  these  states  will  be  described 
in  the  next  report. 


D.  The  Y  State 

The  MC-VCI  energy  of  this  state  is  about  0.25eV  lower  in  energy  than  the 
^A  state.  But  the  natural  orbital  analysis  of  the  7I  state  shows  that  state 
is  of  Fe+2(f2  ionicity  with  the  valence  electrons  of  Fe  in  a  sd5  7SZ  coupling. 

The  s  orbital  is  highly  polarized  which  will  significantly  reduce  the  correlation 
energy  between  the  s  and  d  electrons.  *The  and  Y  state  energies  must  be 
corrected  by  the  differential  correlation  energy  between  the  d6  and  d5  configuration. 
Calculations  are  now  in  progress  to  determine  the  magnitude  of  this  energy  which 
will  probably  exceed  leV. 
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4.  Spin-Orbit  Interaction  in  OH  and  SH,  X^n 

The  spin-orbit  coupling  constant  for  OH,  OD,  and  OT  have  been  reported 
by  Coxon  et  al.  [1,2,3]  as  a  function  of  internuclear  distance.  The  curve 
was  obtained  by  inverting  data  from  the  first  10  vibrational  levels.  The  curve 

O 

has  a  minimum  around  1.5A  as  seen  in  Fig.  1,  and  is  found  to  then  increase 
to  larger  distances.  This  was  quite  unexpected  and  prompted  several  theoretical 
calculations.  An  accurate  all -electron  calculation  by  Langhoff  et  al.  [4]  agreed 
accurately  with  the  experimental  data.  We  had  already  started  a  similar  calculation 
using  the  effective  potentials.  This  was  reported  in  the  37th  Symposium  on 
Molecular  Spectroscopy  held  at  the  Ohio  State  University  in  June  1982. 

The  abstract  of  the  paper  presented  at  that  meeting  is  given  below. 

SPIN-ORBITAL  COUPLING  IN  THE  HYDROXYL  RAO I CAL,  OH 
WALTER  J.  STEVENS  AND  M.  KRAUSS 

Ab  initio  effective  spin-orbit  operators  (AESOP),  derived  from  relativistic 

effective  potentials  based  on  Dirac-Fock  atomic  wavef unctions,  have  been  used  to 

calculate  the  spin-orbit  coupling  in  the  ground  11^  state  of  OH  as  a  function 

2 

of  internuclear  distance.  Spin-orbit  interactions  between  the  n.  state  and  the 
2-4  32 

Z~  and  Z~  states  arising  from  the  0(  P)  +  H(  S)  asymptote  have  been  taken  into 

2  + 

account  as  well  as  interaction  with  the  A  Z  state  which  correlates  with  the 
1  2 

0(  D)  +  H(  S)  asymptote.  The  ground-state  was  represented  by  an  11  configuration 
MC-SCF  wavef unction  while  the  excited  states  were  determined  by  full  valence 
first-order  Cl  calculations  using  the  gound  state  molecular  orbitals.  Results 
near  Rg  are  in  excellent  agreement  with  experimentally  determined  spin-orbit 
splittings.  Spin  recoupling  is  found  to  play  an  important  role  in  the  variation 
of  the  spin-orbit  splitting  with  R.  The  coupling  constant.  A,  does  not  decrease 
monotonically  to  the  oxygen  value  as  previously  thought.  Vibrational  averaging 
of  our  calculated  values  gives  excellent  agreement  with  the  experimentally  determined 
values  of  Coxon  and  Foster  [Can.  J.  Phys.,  to  be  published]  through  v=10. 
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The  vibrational  average  for  the  OH,  OD,  and  OT  molecules  are  shown  in  Fig.  2. 

The  shapes  of  the  calculated  and  experimental  curves  are  in  excellent  agreement 
but  the  experimental  curve  for  OD  is  slightly  lower  than  that  for  OH  which  does 
not  agree  with  the  calculation. 

Using  the  effective  potentials,  it  is  straightforward  to  repeat  these 
calculations  for  SH.  The  vibrationally  averaged  spin-orbit  compling  constants 
for  SH,  SO,  and  ST  are  shown  in  Fig.  3.  The  qualitative  behavior  of  the  OH  and 
SH  spin-orbit  coupling  is  obviously  very  much  the  same.  As  in  all  calculations 
with  the  effective  potentials,  it  is  possible  to  range  widely  over  valence 
iso-electronic  molecules  going  down  a  column  of  the  periodic  table.  The  large 
amount  of  computer  time  required  for  the  all-electron  calculations  precludes 
such  excursions  into  the  heavier  systems. 

The  qualitative  behavior  of  the  spin-orbit  curve  is  given  by  the  spin-uncoupling 
inherent  in  the  base  wavefunction  which  will  dissociate  OH  formally  correctly 
into  fragment  Hartree-Fock  atoms.  The  wavefunction  is  given  in  Fig.  4.  Because 
of  the  open  shell  character  of  both  of  the  atomic  fragments,  the  ba?a  configuration 
has  three  components  including  the  one  with  the  3a4o  orbitals  triplet  coupled. 

This  leads  to  one  of  the  asymptotic  determinants  to  have  a  spin  on  the  hydrogen 
atomic  orbital  opposite  to  that  of  the  other  two.  As  a  result  the  asymptotic 
spin-orbit  coupling  constant  for  the  n  state  is  only  two- thirds  of  the  constant 

O 

of  the  0(  P)  atom.  The  proper  atomic  spin-orbit  constant  is  recovered  by  coupling 
with  the  other  states  that  arise  from  the  0(  P)  +  H(  S)  interaction.  However, 
this  coupling  does  not  become  significant  until  the  energy  separation  between 
the  molecular  states  becomes  comparable  to  the  off-diagonal  spin-orbit  coupling. 

o 

This  occurs  at  distances  of  the  order  of  2.5A  which  is  well  beyond  the  turning 
point  for  v=10.  As  a  result  the  absolute  value  of  the  coupling  constant  will 
decrease  as  v  increases  until  the  molecule  is  excited  close  to  the  dissociation 


limit. 
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